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We consider frequency dependence of the neutron scattering amplitude from a two-dimensional 
quantum antiferromagnet. It is well known that the long range order disappears at any finite 
temperature and hence the elastic neutron scattering Bragg peak is transformed to the quasielastic 
neutron scattering spectrum oc du/uj. We show that the widely known formula for the spectrum 
of an isotropic antiferromagnet derived by Auerbach and Arovas^ should be supplemented by a 
logarithmic term that changes the integrated intensity by two times. A similar formula for an 
easy-plane magnet is very much different because of the Berezinsky-Kosterlitz-Thouless physics. An 
external uniform magnetic field switches smoothly the isotropic magnet to the easy-plane magnet. 
We demonstrate that the quasielastic neutron scattering spectrum in the crossover regime combines 
properties of both limiting cases. We also consider a quantum antiferromagnet close to the 0(3) 
quantum critical point and show that in an external uniform magnetic field the intensity of elastic 
(quasielastic) neutron scattering peak depends linearly and significantly on the applied field. 
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I. INTRODUCTION 

The two-dimensional antiferromagnets have been stud- 
ied thoroughly in numerous works. Nevertheless, sur- 
prisingly, some basic questions related to the elastic and 
quasielastic neutron scattering remain unresolved. It is 
widely known that the (quasi-)two-dimensionality of lay- 
ered systems yields logarithmic corrections to the sub- 
lattice magnetization, amplitude of magnon scattering^ 
and spin correlation functionsi^^. Both, magnetic field 
or easy-plane magnetic anisotropy weaken the above- 
mentioned divergent contributions and simultaneously 
yield a Kostcrlitz-Thoulcss behavior at low temperatures. 
Studying neutron scattering amplitude in the presence of 
these factors is an important problem, which is also rel- 
evant from the experimental point of view. 

Our interest to this problem has been stimulated by re- 
cent experimental studies^ of quasielastic neutron scat- 
tering from undcrdoped YBa2Cu3 6.45- This sample cor- 
responds to the hole doping level about 8.5%. The paper- 
reports a weak quasielastic neutron scattering and hence 
indicates a quasistatic magnetic ordering. This demon- 
strates proximity to a quantum critical point (QCP) be- 
tween magnetically disordered and magnetically ordered 
states. The paper— reports that the quasielastic neutron 
scattering intensity depends substantially on the applied 
uniform magnetic field. Earlier a similar effect of en- 
hancement of the quasielastic scattering in magnetic field 
was observed in Lai.gSro.iCuO^ Ref. [|| We believe that 
the effect in both compounds is of the same origin. How- 
ever, a theoretical analysis in La2_ 2; Sr 2 ;Cu04 is much 
harder due to a large degree of intrinsic disorder related 
to random Sr positions. The observed magnetism is two- 



dimensional (2D), there are no indications for a correla- 
tion in the direction perpendicular to Cu02 planes. Fur- 
thermore, it is incommensurate, mainly due to the con- 
tribution of charge degrees of freedom, which certainly 
complicates a theoretical analysis of magnetic properties, 
see Ref. [l(| Moreover, both compounds, YBa2Cu30g.45 
and Lai.gSro.iCu04 are superconducting at low enough 
temperatures. 

In the present paper we simplify the problem and con- 
sider two dimensional quantum antiferromagnets with or- 
dered ground state described by the Heisenberg model. 
The presented analysis can be also relevant for the lay- 
ered perovskite compounds, in particular K2MnF4 and 
K 2 CuF 4 (Ref. i). We are interested in temperatures and 
external magnetic fields that arc much smaller than the 
Heisenberg exchange interaction J. Therefore, it is quite 
natural to apply techniques of the nonlinear a- model that 
significantly simplifies calculations. It is known for long 
time that the model describes the low-energy properties 
of the collinear quantum magnets^. 

A quantum magnet in a ground state spontaneously vi- 
olates the continuous symmetry of the Hamiltonian and 
this leads to the spontaneous magnetization and, due to 
the Goldstone theorem, to gaplcss magnons. Neutron 
magnetic scattering from the static magnetization leads 
to the elastic Bragg peak. Due to the gapless excitation 
spectrum and due to the low dimensionality, any finite 
temperature destroys the long range order and hence de- 
stroys the elastic Bragg peak. The peak is transformed 
to the narrow quasielastic spectrum. The first problem 
that we address in the present paper is the derivation of 
this quasielastic spectrum. 

The second problem that we address in the present pa- 
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per is how the intensity of elastic (quasielastic) neutron 
scattering depends on the applied uniform magnetic field 
if the system is in the vicinity of a QCP. The obtained 
result is generic and is independent of the specific mech- 
anism of the quantum phase transition. We believe that 
it explains the data obtained in Refs. 0IM 

The structure of the paper is the following: In Sec- 
tion II we formulate problems and present answers for 
quasielastic neutron scattering spectra at nonzero tem- 
perature. In Section III we discuss the dependence of 
clastic neutron scattering on applied magnetic field and 
compare theory with experimental data. Section IV 
presents simple heuristic derivations of quasielastic neu- 
tron scattering spectra at nonzero temperature. A rigor- 
ous derivation of quasielastic spectra based on the Renor- 
malization Group (RG) analysis is given in Section V. 
Section VI presents our conclusions. 



II. QUESTIONS AND ANSWERS 

In the present Section we only formulate problems 
and present our results for neutron scattering spec- 
tra/probabilities in various cases. Derivations of these 
results are presented in sections IV and V. 

The Lagrangian of the nonlinear 0(3) a- model describ- 
ing a 2D isotropic quantum antifcrromagnet reads^ 
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where x_L is the magnetic susceptibility and p s is the 
spin stiffness. The unit vector, n ~ (n x ,n y ,n z ), n 2 = 1, 
describes staggered spins. Here we consider "quantum 
renormalized" model with x± and p s being renormalized 
(observed) parameters of the ground state. We use the 
choice of normalization (cutoff), such that (rt z )y =0 = 1. 
The ground-state sublattice magnetization p is consid- 
ered as an independent parameter of the theory^. Also, 
we consider real time since our aim is the real excitation 
spectrum. Throughout the paper we set the Planck's 
constant and the Boltzmann's constant equal to unity, 
H = fee = 1- Interaction with the incident neutron we 
take in the following simplified form that is sufficient for 
our purposes 



(2) 



where a is the Pauli matrix describing the neutron spin, 
and Vv is the wave function of the neutron. Since ft is the 
staggered magnetization, Eq.([2]) implies that the momen- 
tum transfer is shifted by the AF vector. We also intro- 
duce external magnetic field by changing ft — > ft— 

At zero temperature the 0(3) rotational symmetry is 
spontaneously broken and the ground state magnetiza- 
tion is Hq = (0,0, 1). Quantum fluctuations around this 



ground state are, 
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Here e\ (A = x, y) is the magnon polarization, the 
unit vector along the corresponding direction in the spin 
space; a kx is creation operator of the corresponding 
magnon; and V is the area of the system. 

Application of the Fermi golden rule to the Hamilto- 
nian ([2]) with account of ((3]) and the averaging over neu- 
tron polarizations gives the elastic and inelastic scatter- 
ing probabilities with momentum transfer q and energy 
transfer u) 



W(u, q) = V [8{io)5 qfi + 2A 2 q 5{u - u q )] 



(4) 



Here 5 qt k is the Kronecker symbol. This gives the well 
known result for the q-intcgrated scattering probability 
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The subscript shows that this is the zero temperature 
case. 

The 0(2) nonlinear a model describes an easy-plane 
AF. The Lagrangian is the same ([T]), but the spin vector 
is two dimensional, n = (n y ,n z ). Obviously, in this case 
the elastic scattering probability is the same, and the 
inelastic one is twice smaller because there is only one 
magnon. 
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It is well known that any small but finite temperature 
destroys the magnetic ordering, so the 5-function peaks 
in both Eqs. ([5]) and ([6]) disappear, they are replaced 
by broad quasielastic peaks. In the 0(3) case a nonzero 
temperature leads to a finite correlation length £ and to 
the corresponding "spin- wave gap"— 
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We have set the lattice spacing equal to unity and intro- 
duced the dimensionlcss temperature 
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In the 0(3) case our result for the q- integrated neutron 
scattering probability at nonzero temperature reads 
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The probability is sketched in FigJT] by the black solid 
line. The nontrivial point is that the quasielastic peak at 




FIG. 1: (Color online) A schematic sketch of the momentum 
integrated neutron scattering probability Wt{w) at nonzero 
temperature. The probability for an isotropic magnet (0(3) 
(j-model) is shown by the black solid line. The probability 
for an easy-plane magnet (0(2) a- model) is shown by the red 
dashed line. 

M ^ T differs from the previously known formula^ by 
the logarithmic term -rduj/\u\ ln(T/|w|)/ln(T/A). The 
term is important to satisfy the spectral sum rule: the 
total quasielastic intensity at nonzero temperature must 
be equal to the elastic intensity at zero temperature 
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This sum rule must be valid with relative accuracy up to 
terms of the order of ~ r <C 1 . Without account of the 
logarithmic term the sum rule is violated. It is easy to see 
that the term reduces the total quasielastic intensity by 
two times. A simple heuristic derivation of the spectrum 
([9]) is given in section IV, and a rigorous derivation based 
on the RG analysis is given in section V. 

An easy-plane quantum magnet is described by the 
0(2) nonlinear cr-model. The long range magnetic order 
is also destroyed by a nonzero temperature, but in this 
case there is no a correlation length and no a "spin-wave 
gap" . Because of the Berezinsky-Kosterlitz-Thoulcss 
physics all correlators decay as powers of distance^. Our 
result for the nonzero temperature q-integrated neutron 
scattering probability in this case reads 
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The probability is sketched in FigJT] by the red dashed 
line. The energy scale A plays a role in this case as 
well. Certainly it is not a gap any more, this is just 



the relevant energy scale. At |w| 3> A the intensity (fTTj) 
is twice smaller than ([9]). This is because there is only 
one magnon in the 0(2) case against two magnons in 
the 0(3) case. On the other hand at |w| < A the 0(2) 
intensity is much larger than the 0(3) one. The spectral 
sum rule (|10|) is satisfied in this case as well because 
the elastic intensity is the same in both 0(2) and 0(3) 
cases. A heuristic derivation of the spectrum is given 
in section IV, and a rigorous RG derivation is given in 
section V. 

A uniform magnetic field, applied to an isotropic an- 
tifcrromagnct, orients sublattice moments perpendicular 
to the field and therefore acts as a kind of easy plane 
anisotropy. This means that the magnetic field smoothly 
transfers the 0(3) magnet to the 0(2) magnet. If temper- 
ature is zero then the effect of the magnetic field acting 
on isotropic antifcrromagnct is very simple and well un- 
derstood. We consider here a nonzero temperature. For 
convenience we rescale the magnetic field, g^sB — > B, 
where g is the g- factor, and /is is Bohr magneton. If the 
magnetic field is larger than temperature, B > T, then 
obviously the quasielastic peak is described by the 0(2) 
formula (jlip . If the magnetic field is very small, B <C A, 
then practically the quasielastic peak is described by the 
0(3) formula ©. The only nontrivial case is the case of 
the intermediate magnetic field, A <C B <C T. Our result 
for the q-intcgratcd spectrum in this case reads 
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The probability (fT2|) is sketched in Figj2j Naturally, the 
spectrum (fTSj) satisfies the spectral sum rule (fTO|) . A 
heuristic derivation of the spectrum (|12[) is given in sec- 
tion IV, and a rigorous RG derivation is given in section 
V. 



III. DEPENDENCE OF 
ELASTIC/QUASIELASTIC NEUTRON 
SCATTERING INTENSITY ON MAGNETIC 
FIELD IN A VICINITY OF A MAGNETIC 0(3) 
QCP 

In this section we discuss enhancement of the staggered 
magnetization by applied external magnetic field. This 
effect was first pointed out for the square lattice Heisen- 
berg antiferromagnet in Refill In the present paper we 
stress that the effect is strongly enhanced in the vicinity 
of a magnetic QCP. To discuss this effect it is sufficient to 
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FIG. 2: A schematic sketch of the momentum integrated 
neutron scattering probability Wt(w) for isotropic antiferro- 
magnet (0(3) a-model) in the external magnetic field in the 
regime A « fi « T. 



consider the zero temperature case. At a finite temper- 
ature the integrated quasielastic intensity equals to the 
elastic intensity at zero temperature. The magnetic field 
changes the staggered magnetization, 



(14) 
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For further numerical estimates we refer to the cxper- 
imenti performed with YBa2Cu3 06.45- The experimen- 
tal dependence of the quasielastic neutron scattering in- 
tensity on magnetic field is shown in Figj3] Note that 
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FIG. 3: (Color online) YBa2Cu306.4s: The experimental de- 
pendence^ of the elastic peak intensity on the external mag- 
netic field at T = 2 K. The intensity is measured at the in- 
commensurate position of the neutron scattering peak. The 
line is the result of a linear fit to the data points. 



We remind that our normalization is (?t. z )t=o.b=o 
1 . In the presence of external magnetic field the in- 
plane magnon remains gapless (easy plane), ujk = ck, 
while the out-of-plane magnon becomes gapped^, ujk = 
y/c 2 k 2 + B 2 . This results in the enhancement of the stag- 
gered magnetization, 
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Clearly this formula is valid only if 6(n z ) <C 1. There 
are two important points to note about Ea.([j"5"]): (i) the 
dependence on the magnetic field is linear, (ii) the de- 
pendence is enhanced approaching a QCP where the spin 
stiffness p s vanishes (The Eq. (| 15[) is however inapplica- 
ble in the near vicinity of the QCP where the criterion 
5{n z ) <C 1 is violated), see e.g. Ref. ITil. 

To estimate magnitude of the effect (fl~5j) we refer to 
the simple models, two coupled Hcisenberg planes with 
spin 1/2, the position of the QCP is g c = J±/J = 2.525. 
In the vicinity of the QCP the observed staggered mag- 
netization p and the spin stiffness scale as p cx (g c — , 
p s oc (g c — g) u , where the critical indexes are related, 
2/3 = (1 + T])u, see e.g. Refs. HE The index 77 w 0.03 is 
very small, so practically v « 20. The spin stiffness at 
J± = is p s0 w 0.2 J and the staggered magnetic moment 
p 2D w 0.6pB, see Ref. [H. Hence 



the compound is a superconductor with critical super- 
conducting temperature T c = 35K. This is the second 
order superconductor, the magnetic field of several Tesla 
pretty much penetrates in the bulk, and therefore we 
disregard superconductivity in the analysis. The com- 
pound is very close to the magnetic QCP, the value of 
the static "staggered" moment is p ~ 0.05^b, see Ref. 0- 
Therefore, according to Eq. (fT6]) the spin stiffness is about 
p s ~2x 10~ 3 J. With the magnetic field B = 15 Tesla 
and with the characteristic cupratc exchange integral 
J = 1400A", Eq. (TT5|) gives the following field induced 
variation of the magnetization 



Sp 



\B\ 

&np s 



0.3 



(17) 



This value is very close to the datai presented in FigJ3] 
(Intensity is quadratic in magnetic moment, therefore 
SI 1 1 = 2Sp/p). There is no doubt that the magnetic 
QCP observed in YBa2Cu30 w is of a more complex na- 
ture than just a simple 0(3) QCP, see a discussion in 
Ref. ITol. Nevertheless, it is clear that generically Eq. (fT"5|) 
must be valid because it is based only on two physical 
inputs: gapless spin-wave spectrum and reduction of the 
effective spin stiffness in the vicinity of the QCP. 

The La based compound La2- a; Sr a ;Cu04 is very 
strongly disordered due to random positions of Sr ions. 
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FIG. 4: (Color online) Lai.gSro.iCuCU: The experimental 
dependence^ of the static magnetic moment on the external 
magnetic field at T — 2 K. The magnetic moment is mea- 
sured at the incommensurate position of the neutron scatter- 
ing peak. The data points shown by red circles are obtained 
from Fig.3a of Ref. [8|, and the data points shown by blue 
triangles are obtained from Fig.3b of Ref. @. 



In addition the compound has a nonmonotonous de- 
pendence of the static magnetic moment on dopingii 
(anomaly at x « 0.12). These two reasons complicate 
a comparison between theory and experiment. The mag- 
netic field dependence of the static magnetic moment in 
Lai.gSro.iCuO.4 has been measured by neutron diffrac- 
tion.— The experimental magnetic moment versus mag- 
netic field is plotted in Fig HI The magnetic moment is 
normalized to its value at zero field. FigfJ] displays the 
quick growth from zero to 5 Tesla and then a slow lin- 
ear dependence. A possible explanation of these data 
is based on the idea of magnetism induced around su- 
perconducting vortex coreJ^ An alternative explanation 
based on the approach developed in the present paper is 
that the steep dependence up to 5 Tesla is due to glassy 
effects related to random positions of Sr ions, and the 
slow dependence for 5 < B < 15 is due to suppression of 
quantum fluctuations described by Eq. (fT5]l . Within this 
framework the effect observed in both YBa2Cu3 06.45 and 
Lai.gSro.iCu04 has the same physical origin. 



IV. HEURISTIC DERIVATION OF INELASTIC 
NEUTRON SCATTERING SPECTRA AT 
NONZERO TEMPERATURE 

We start from derivation of Eq.© that describes 
quasielastic scattering from the isotropic antiferromag- 
net (the 0(3) case). Let us consider the energy trans- 
fer uj much larger than the "spin wave gap" and some- 
what smaller than temperature, A <C < T. Hence 
the momentum transfer is q ~ ui/c. The correspond- 
ing time/space scales, 8t ~ 1/w, 6r ~ l/q, are much 
smaller than the time/space scales of slow thermal fluctu- 



ations that destroy the long-range AF order, Stf ~ 1/A, 
Srf ~ £. Thus, the neutron makes a snapshot of the 
slowly fluctuating antiferromagnct, the AF order is un- 
touched at the relevant scale. Hence, the formula (U]) for 
the scattering probability is almost correct. One needs 
only to "erase" the elastic peak, and in the inelastic part 
to account for thermally excited magnons with the stan- 
dard Bosc occupation numbers 
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e <WT - l 
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Hence, ((4]) is transformed to 



W(u, q) = 2VA 2 q [(n q + l)S(u> - u q ) + n q< 5(w + w,)](19) 

As usually the first term comes from the magnon emission 
and the second term comes from the magnon absorption. 
Integrating (|19[) over q and expanding at \ui\ <C T one 
finds 



Wt(uj)cIuj —> tV 



du) 



(20) 



To propagate the spectrum down to ui ~ A we note that 
the power of uj is fixed by dimension. However, dimension 
does not forbid an additional logarithmic contribution, 
Wt -t Wt [1 + ^m(T/|u;|)], where fi <C 1 is some small 
coefficient. The logarithmic term must be negligible at 
T > \uj\ > A, where Eqs.([T9 )l .([2"0" ]) have been derived, 
however it can be 100% important at \cj\ ~ A. Impos- 
ing the spectral sum rule condition (|10[) one immediately 
finds that fi rts — l/ln(T/A) w — r. This concludes the 
heuristic derivation of the spectrum ([9]). Notice that the 
spectral sum rule uniquely determines the spectrum as 
soon as one takes the first power of logarithm. However, 
a priory, one cannot exclude higher powers, for example, 
W T -> W T [l + Mln(T/M) + ^ln 2 (T/|u;|)] . Therefore, 
the presented simple heuristic derivation certainly cannot 
substitute a rigorous derivation. The rigorous derivation 
of © based on RG is given in Section V. 

In the 0(2) case (easy-plane magnet) the logic of 
derivation is very similar. In this case Eqs. (|19[) and 
(|2TJ1) must be also valid in the regime A <C w < T. The 
only difference is that the prcfactor in these Eqs. must 
be two times smaller because in the 0(2) model there 
is only one magnon instead of two magnons in the 0(3) 
model. However, propagation of these Eqs. to smaller uj 
is different because of the Berezinsky-Kosterlitz-Thouless 
physics in the 0(2) case. In this case there is no a corre- 
lation length, all correlators decay as powers of appropri- 
ate variables^. Therefore, the only possible correction 
to the spectrum (|20|) is an additional prefactor (T/ui) p , 
where p <C 1. The quasielastic spectrum still has to sat- 
isfy the spectral sum rule (jlOp . From here one immedi- 
ately finds that p = — r and this proves the validity of 
the spectrum (|11[) . An RG derivation of the spectrum 
(TTT1) is presented in the Section V. 

Eq. (fT2")) describes the quasielastic scattering from the 
isotropic antifcrromagnet in the external magnetic field, 



6 



A <C B <C T. Obviously, at a frequency higher than the 
magnetic field, B < \uj\ < T, the field is not important 
and hence the spectrum (fT2|) must coincide with ([9]). At 
|cj| < B the gapped magnon is switched off, the prob- 
lem is effectively reduced to the 0(2) case, and hence the 
shape of the spectrum must be the same as that in (fTT|) . 
However, one has to take care about thermal fluctuations 
with typical frequencies between B and T . These fluctu- 
ations reduce value of the magnetization that comes to 
the effective 0(2) cr-model regime, 



(n z ) 2 = 1 -> (n z ) 2 = l-(n 2 



(21) 



Neutron scattering is proportional to the magnetization 
squared, therefore, (f2"Tj) presents another explanation of 
Eq.©. Because of the reduction of number of active 
magnons, 2 — > 1, the intensity at w = B — must 
jump down by factor 2 compared to the intensity at 
cu = B + 0. Further evolution down in frequency is 
governed by the 0(2) cr-model. However, we have de- 
fined the effective 0(2) cr-model with (n z ) = 1 at the 
upper cutoff which is B in this case instead of T in (fTTI) . 
Therefore to apply one has to renormalize magne- 
tization back to unity, (n z ) — > 1, and this implies that 
p s — > p* s = p s [l — t In (5)1 ■ This explains why r defined 
in (|8) is replaced by t* given by (fL3")l . Notice that the 
spectrum (fT2]) satisfies the spectral sum rule (| 10|) . Rig- 
orous RG derivation of this spectrum is presented in the 
Section V. 



RENORMALIZATION GROUP (RG) 
ANALYSIS 



To perform calculations of thermodynamic and dy- 
namic properties we consider contribution of different 
momenta in small steps. To account for thermal fluc- 
tuations, we put ultraviolet cutoff A on the theory 
(which is temperature T/c — > T in the finite-temperature 
renormalized-classical regime) and follow the standard 
procedure 2 *^. Supposing that magnetic field acts along 
the x-axis and the order parameter is along the z-axis, we 
represent the vector n = (% x , tt v , yl — tt 2 ) and introduce 
the two-component vector n =(ftx,'Ky). In terms of the 
two-component field 7r the Lagrangian becomes (now we 
use the imaginary time t, and set the spin-wave velocity 
equal to unity, c = 1) 

C = 2L J dtd 2 x[(diT) 2 + B 2 ir 2 , + hn 2 

+ {nd7r) 2 + ^tt 4 - 4iBTT 2 d t TT c - B 2 



-(n 2 + 1 -^)} 
Ps 2 



(22) 



where d = (V,9 t ), and we have introduced the source 
field h for the longitudinal (n z ) spin component. The 
terms in the first line of Eq. (f2"2"f correspond to the 



Lagrangian of the noninteracting spin waves, the other 
terms are responsible for the spin-wave interaction (the 
last term arises from the condition n 2 = 1). 

Let consider the renormalization of parameters of the 
spin wave spectrum by the interaction. To first order in 
1 1 p s we have 
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After performing decouplings (|23|) we obtain up to a con- 
stant contribution 



^ / ' dtd 2 x\Y^. (d^fil + in 2 )) (24) 

/ — ^i—x.y 



Ps 
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+B 2 ir 2 (l - (<)) + hir + 2Bn 2 d t -Ky + 0(ir 4 



where we have taken into account ((diTi) 2 ) — l/p s = 
—m 2 (ir 2 ); m 2 = B 2 + h, m 2 = h arc masses of ir x 
and TT y fields respectively. Rescaling the field tt —> 
(1 + (7r|)) _1 / 2 7r to restore the coefficient at (dn) 2 we ob- 
tain 



C 



% jdtd 2 x [(dn) 2 + B 2 7 rUl-2(n 2 x )) 
+2BTT 2 d t n y + 0(tt 4 ) 



(25) 



where x = y and y = x. The corresponding renormalizcd 
parameters are 



B 2 (l~2(n 2 x )) 



h(l 



1 



<*?» 



(26) 



From the performed rescaling we have in addition 

p*Z 2 =p 8 (l + (n 2 )) (27) 

where Z^ is the rescaling factor of the field . Requiring 
similar to-^ h R pf i = Zi(hp S i), we obtain 



Therefore 



Pi 



Zi = 1 



(28) 



p.0- -{*§)) 
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The renormalization factor for the n z field can be deter- 
mined from the sublattice magnetization 



(»«) 



1 



1 



(30) 
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which implies Z x = Z y = Z z . Similarly to the isotropic Representing I = ln(A//i) we obtain 



case*^, the obtained expressions for the renormalized pa- 



rameters (|29f can be used to formulate the rcnormaliza- 
tion group transformation. Below we consider the renor- 
malized classical regime, T 3> max(A, B). 

As it is argued in Refs. HQ it is sufficient in the renor- 
malized classical regime to account only for classical (zero 
Matsubara frequency) contributions to (7rf). Consider- 
ing the contribution of infinitely thin layer in momentum 
space, we obtain (r^ = T/(2np t s )) 



dl 

dTy_ 
dl 

d\nZ 2 
dl 



1 + e 2l B 2 ' 
1 + e 2l B 2 ' 



(31) 



where I = ln(A//x) is the scaling parameter. The flow 
of the magnetic field can be neglected in the one-loop 
approximation. 

For B — Eqs. (|3"Tj) are equivalent to standard one- 
loop RG equations for 2D antiferromagnet^iiS. Eqs. (f3Tj) 
are similar to those for the easy-plane case studied in 
Ref. H3, although differ from those due to presence of 
two spin stiffnesses instead of one. The other way of 
writing first two Eqs. (|3ip is 



dhiT x 

dl 
d In T y 

dl 



Ty, 



1 + e 2l B 2 



(32) 



Together with the last Eq. (|31j) this implies that 
^ 2 (0/[ r *(0 r 2/(0] = I/ 7 " 2 i s the invariant of the flow. 

At e 2l B 2 a (l) < 1 [0(3) regime] we have t x (1) w t v (1) 
such that 



di~ x , y 2 

dl ~ Tx > v 

and T XyV {l) I 'Z(l) = r, which yields 

Tx (l) « r s (i)«r/(l-rO; 
Z{1) » 1/(1 -rZ) 



(33) 



(34) 



At e 2l B 2 (l) ~ 1 the crossover from 0(3) to 0(2) regime 
occurs. At the crossover scale Z* = ln(A/_B) we obtain 



tZ* with 

Z* w 1/[1 -rln(A/B)] 



(35) 



For e 2l B 2 a (l) > 1 [0(2) regime] the flow of r c (Z) is sup- 
pressed and we have r y (l) = T^°=const, 



eZlnZ 2 



din 7% 



2tt 



where r 



dl dl 
— t* up to logarithmic accuracy, such that 

r x (Z) « r* exp[r* (Z-Z*)] 
Z 2 (l) « Z* 2 cxp[r*(Z-Z*) 



(36) 



T"x(M) 

Z 2 ( M ) 



rZ*(fl/At) 1/ln(i,/A) 



(37) 



where Z* is given by Eq. ([33)) and to one-loop order A = 
Aexp(— 1/t) (Note that the two- loop corrections yield A 
given by Eq. Q , see Ref H|) . This power-law depen- 
dence corresponds to the Berezinsky-Kosterlitz-Thouless 
behavior of correlation functions (sec below). 

Now we discuss the behavior of correlation functions 
important for experimental data. From our scaling anal- 
ysis we can obtain the transverse Green functions 



(X,K-q}}^ = 



1/T 



dt(T[nUt)nU0)]) ( 



(38) 



(«k! q »< 



2-KTi, 



iq 



^/q 2 + to 2 , m x = B, m y 



0; Tiq , Z^q 



ilYV 



where u>i, 

obtained from the result of the solution of Eqs. (|3"Tj) by 
substitution Zj = ln(T/aj.; q ). Therefore, we obtain the 
results for the functions at |q| ^> B in the renormalized 
classical regime 



(Kl»-q»*- = (K\n%))i 



2tTT* 



1 



Ifl-rlnA 



TZ 2 q 2 
1 



(39) 



This result agrees with the result for the correlation func- 
tion derived for the 0(3) model by Chakravarty, Halperin 
and Nelson^. On the other hand, in the scaling limit 
Z — > oo we obtain from Eq. ([3T1) 



1 - r In ■ 



A 



l/ln(B/A) 



(40) 



Z 2 ( M ) V" "B 
which implies the results for the functions at |q| <CB 



i 

Ps 



1 - t In 




V 



(41) 

with a = l/ln(_B/A). In this regime we obtain the 
power law corrections to the Green functions, which are 
in agreement with the Berezinsky-Kosterlitz-Thouless so- 
lution of the 0(2) model. Note, that both results, ([31?]) 
and (1411) fulfill the sumrulc 



[(KK)) U + KK»J = 1 ( 42 ) 



which is actually required for the total spectral weight 
including the longitudinal Green function. This shows 
that the contribution of the longitudinal Green function 
is contained implicitly in the calculated transverse Green 
functions, which seems to be the general property of the 
1/5 expansion^!. 
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For the structure factor, given by imaginary parts of 
the Green functions at the real axis, we obtain 



7T ^ — ^■i—x,y 



(43) 



' ?q 



[S(tJ - CJiq) + S(uj + CJ lq )] 



where n u is the Bose function. Integration over q yields 

~~ ,,,2 



/ jH<K"-^)+%+^)] 

— — - > v(oj — mj) — 7t 



(44) 



where Ti W = t;(7 = ln(T/|o;|)) and similar for Z u . 

In the 0(3) regime \w\ > B where T xu] /Zl = T yu /Z% 
t[1 — Tln(A/|w|)], we obtain 



t ( A 
W(e> ) = T - l-rln n 



(45) 



which agrees with the result ©• In the 0(2) regime 
Qw\<B) 



zi 

Zl 



const 








r(l-rl»A) 









(46) 



Bragg peak observed in neutron scattering. At an ar- 
bitrary small but finite temperature the long range or- 
der is destroyed and hence the elastic Bragg peak is 
transformed to the quasielastic spectrum. We derive 
these spectra. The derived quasielastic spectrum for 
an isotropic antifcrromagnet is given by Eq.Q and the 
corresponding physics is determined by the finite cor- 
relation length. The derived spectrum §§§ differs from 
the previously known answer— by a logarithmic term 
that changes the integrated intensity by two times. The 
derived quasielastic spectrum for an easy plane anti- 
ferromagnet is given by Eq. (|lip and the correspond- 
ing physics is determined by the Berezinsky-Kosterlitz- 
Thoulcss power law. An external uniform magnetic field 
drives a crossover from the isotropic antiferromagnet to 
the easy plane one. Hence the quasielastic spectrum 
evolves with the magnetic field. This evolution is de- 
scribed by Ea.(|T2]). 

An external uniform magnetic field, applied to an 
isotropic quantum antifcrromagnet, suppresses quantum 
fluctuations. As a result the integrated intensity of clas- 
tic (quasielastic) neutron scattering grows linearly with 
applied external magnetic field, see Eq. (fl"5")) . The de- 
pendence is significantly enhanced in the vicinity of the 
quantum critical point that separates magnetically or- 
dered and magnetically disordered states, an estimate is 
given in Eq. pT)) . We believe that the observed enhance- 
ment of the quasielastic scattering in magnetic field in 
YBa 2 Cu3 06.45, Ref. 0, and in Lai.gSr .iCuO4, Ref . H, is 
due to this mechanism. Extension of the present analysis 
to include the magnetic incommensurability and contri- 
bution of conduction electrons is of a great interest. 



and 



in agreement with (|12l) . 



(47) 



VI. CONCLUSION 

We have considered two dimensional quantum antifer- 
romagncts. At zero temperature the antiferromagnets 
possess the long range order that results in the clastic 
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